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ABSTRACT 


Light falling on an engineering surface will be reflected partly 
in a preferred direction (specularly), partly in a diffuse manner. The 
scattering of the light energy is assumed to follow Lambert's cosine 
rule. The force differential acting on an elemental surface is set up 
for such cases of mixed reflection. A smoothness coefficient (g 1) 
accounts for that portion of the reflected energy that is radiated 
specularly. A loss coefficient, also § 1, gives the fraction of the 
incident energy that enters into the reflected light. 

Applications are made to the sphere and its motion under the impact 
of light, and to the circular cone with fully illuminated curved surface, 
some results being added when a shadow region exists on it. 
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TECHNICAL MEMORANDUM X-53617 


ON MIXED REFLECTION OF SUNLIGHT 
SUMMARY 


If an elemental surface is exposed to sunlight the total light force 
on it is the compound effect of three intensities of which that of inci- 
dent light can be taken as the basic measure. Only a portion of it will 
be reflected (described by a loss coefficient, k). This portion in turn 
will usually be radiated back in part specularly, in part diffusely. A 
smoothness coefficient, q, identifies the fraction of it that follows the 
familiar mirror reflection law. Lambert's cosine rule is assumed for the 
energy flux distribution in diffusion. On this basis an expression for 
the differential of the light force is derived which must be integrated 
over the illuminated part of the surface. This is more or less easily 
done if the surface equation is known in analytic terms and if the loss 
and smoothness coefficients can be taken as constants. 

An application to the sphere shows that, unless diffuse reflection 
is present, the force on it is that of incidence alone no matter what 
portion of the energy is reradiated specularly. A non-homogeneous sphere 
will, aside from following a straight translation course, in general 
carry out a swinging motion about an axis through the mass center, except 
in singular circumstances when the initial torque is zero so that a 
rotatory motion cannot develop. 

The force on a fully illuminated circular cone can be larger or 
smaller than the incident force alone, depending mainly on whether the 
opening angle of the cone is large or small. The force moment will 
tend to align a homogeneous cone's axis with the light direction, pro- 
vided that the product q K is smaller than 0.6. If it is larger, suf- 
ficiently slender cones are unstable. (The Apollo capsule is likely to 
be stable if the conical part is fully illuminated.) 

A numerical computation was carried out for loss-free reflection 
(K = 1), full illumination, and the opening angle of the Apollo capsule. 
If the axis inclination towards light direction is increased* the force's 
inclination increases also, and does so at a faster pace with smoother 
surfaces. With these, it also becomes larger in magnitude. At a fixed 
attitude, the force depends on the surface quality, growing stronger 
and less inclined as diffusion grows more preponderant. 

In a concluding section approaches other than Lambert's rule to 
account for diffuse reflection are briefly- discussed. 


I. INTRODUCTION 


The energy of reflected light is smaller as a rule than that of the 
incoming light. Part of the latter may go into different forms (heat; 
electrical energy), another part may be transmitted through a trans- 
parent surface. All losses of radiant energy will be consolidated into 
one fractional loss coefficient, K, so that 


I 


r 


K I. 


x 


( 1 ) 


is the fraction of the incident intensity that is preserved in the 
reflected light (0 S K S 1). 

Except in ideal circumstances the surface is rough relative to 
photon size. Experience shows that a portion only of the reflected 
energy resides in light reemitted specularly, that is* in such a manner 
that the surface normal bisects the angle between the directions of 
incoming and reflected rays. A smoothness coefficient q(0 § q § 1) will 
be used to describe the intensity, I s , of light reflected specularly: 


With q = 1, the reflection is ideal; with q = 0, it is diffuse in its 
entirety. Neither of these limiting values will presumably be found 
with engineering surfaces. 

Diffuse reflection is usually taken as obeying Lambert’s rule 
(which, by the way, stems from 1760). The total solid angle (0 = 2it) 
above an elemental surface is divided into elements dft surrounding the 
surface normal. If one of these elemental solid angles makes the angle $ 
with the surface normal, the flux of the light energy in it is assumed 
proportional to cos -9. It may be written as 


d 2 0 = (A cos 9 dfl) dS 


where dS is the area of the surface element and A is a factor of propor- 
tionality. For the integration one best uses the system of spherical 
coordinates cp, -9, r as depicted in figure 1. In it 


dQ, = sin ■O d-9 dcp. 


2 



Figure 1. Spherical Coordinates 
so that the total flux given out by the element dS becomes 


0 


t 



cp=0 -3=0 


2 

sin ■S cos d-9 dtp = rtAdS. 


( 3 ) 


Note that it is proportional to dS. 

In being reemitted light moves in the general 

£ = sin -d cos cp+j, sin •O sin cp + k cos 


direction (figure 1) 

* 




Vectors are denoted by underlined symbols. For example, jL, j., k are 
the unit vectors in the x,y,z-directions , respectively. 
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it exerts a force whose direction is opposite to £ and whose magnitude 
is given by 


1 

c 


d 2 0, 


where c is the speed of light. The total force acting through diffuse 
reemission on the element dS becomes 




2n jc/2 



cp=0 -3=0 


sin -9 cos -9 cl 6 dcp (i sin -9 cos cp + j. sin -9 sin cp 


+ k cos -9) . 


On integrating over cp the and j-components cancel. Thus, 

n/ 2 

k dS J sin -9 cos 2 9 d-9 = - k dS. 

9-0 


The "diffusion" force is in the direction, -k, of the interior normal of 
the surface element. Its absolute value, 


2nA 


dS, 


can be connected with an energy flux normal to dS of magnitude 


„ 2t0\ , c 2 

0 d ■ T dS 3 V 


according to expression (3). On dividing through by dS one obtains an 
intensity related to f,: 



The unknown factor A is eliminated in this expression. Note that while 
in fact the entire fraction I t = (1 - q) I r of the reflected intensity 
is emitted diffusely, only part of it contributes to the accelerating 
force on dS, since the lateral force components cancel (and thus may be 
said to merely produce stresses in the surface element). This is con- 
sequent on the assumption made in the derivation that Lambert's intens- 
ities give rise to actual forces attacking the element in the directions 
-e. 


II. THE LIGHT FORCE ACTING ON DIFFERENTIAL AND INTEGRAL SURFACES 


The intensity of parallel (sunlight) radiation is numerically iden- 
tical with the pressure of light on a surface placed normal to ray direc- 
tion. For such surfaces, expressions (1), (2), and (4) may also be written 
as 


p r = ^o 

P s = <qp 0 > 

p d = I K(1 - q) p 0 . 


(5) 


The pressure, p Q , of the incident light must be taken from observation. 
In nearby space its value is very small: 


P 


o 



( 6 ) 


This figure is an average, for in fact, it varies slightly with the 
earth's distance from the sun. 

An element, dS, sitting on an extended surface will not in general 
be perpendicular to ray direction. If its normal makes the local angle, 
a, of incidence with that direction, the magnitude of the force on it 
is calculated as the product of the pressure with the projected area, 
dS cos a. With incident parallel light the direction of the force is 
that of the rays which we may designate by the unit vector s_. Thus, the 
force element engendered by the light incoming at dS is given by 


d 2 J\ = p Q dS cos a s_* (7) 
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that is, the force is largest at normal incidence (a = 0). It is zero 
with grazing incidence (a = 90°). This is not strictly correct with 
rough surfaces where a few peaks will still be illuminated. However, 
their projected area is small at least of second order. The effect 
can be disregarded (here as well as with a ^ 90°' 

With specular reflection, the force is in a direction, s_' , opposite 
to that of the outgoing light and again makes the angle a with the 
interior normal unit vector n' . Since s_ and s ' are also unit vectors, 
the diagonal's length in figure 2 is 2 cos a. It follows that 


s ' = 2 n cos o: - s . 



Figure 2. The Direction s_' 


Finally, the pressure in the outgoing light being p s , the elemental 
force produced by it becomes 


d 2 F = qK p cos a dS (2n cos a - s). ( 8 ) 

—so — 

Earlier it was found that, with diffuse reflection, the force points 
into the direction of the interior normal (a = 0°). The third relation 
(5) then yields 


d 2 ^ = 3 K(1 - q) P Q dS n 


(9) 
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Expressions (7), (8), (9) add up to give the total force acting on 
the surface element: 


d 2 F = P Q dS 

This formula simplifies considerably in the limiting cases of no reflec- 
tion at all (k = 0) and of ideal specular reflection (< = q ■ 1), less 
so in the case of purely diffuse reflection (q = 0). It will be noticed 
that even with grazing incidence (cos a = 0), there is still a non- 
vanishing force. It is entirely produced by Lambert reflection, dimin- 
ishes with it and disappears in the ideal case q = 1. One may argue 
that with strongly oblique incidence diffuse reflection is not likely 
to follow the rotational symmetry of the cosine law. This possible 
defect, however, would seem to carry weight with plane surfaces only 
where such incidence will extend over finite areas. 

For application of the differential (10) the surface in question 
is best given a pointwise representation: 


s(l - q/c) cos a + 



+ q cos 


-a) . 


( 10 ) 


X = x (o’, t) , y = y(a, t), z = z(a,x). 


( 11 ) 


The intervals in which the parameters a and t move are determined by 
the geometrical shape of the surface. On it, the curves o - const. 
t = const, are coordinate lines (not necessarily orthogonal to each 
other). If the surface is composite, several of such sets (11) will 
have to be used. 

Once expressions (11) are known the components in the (x,y, z) -system 
of the surface normal at any point (a, t) are determinable: 


z a 

1 

Z 0 X 0 

^ 1 

X 

Q 

Q 

z t 

• "2 ■ * N 

Z T X T 

. n 3 * ± N 

X T Z T 


with 


N 


z o 

y-v z t 


z 0 x cr 
Z T X T 


y 

o J o 

x 2 
T T 


( 12 ) 


(13) 
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The indices here signify partial differentiation (y Q = 5y/5a, etc.)* 
From geometrical inspection, the signs in formulas (12) can be chosen 
such that the vector n = n^ + n^ 4- n 3 k points toward the interior, 
that is, toward the region shielded from light rays by the elemental 
surface itself. 

The surface differential dS is given by 


dS = N da dx . 


(14) 


In integrating, da and dx must be defined as positive increments, because 
N and dS are positive. 

All geometric quantities appearing in the force differential (10) 
are known from expressions (12), (13), (14), including the cosine of the 
local angle of incidence, as it is the scalar product 


cos a = n * s . 


(15) 


The light direction s^ is of course a given unit vector. 

Shadow boundaries will often exist on the surface; the integration 
then may use only those parts of the a- and x-intervals that describe the 
illuminated regions. More detail on this can be found in Aero-As trodynamics 
Research Review No. 4. 

Calculation of torques requires the determination of the center of 
pressure for which a general method is set forth in the same paper.* 

By way of examples, the differential (10) will be integrated for 
the sphere and the circular cone. The surfaces are considered non-trans- 
parent to avoid the complications arising through backside reflection. 
Furthermore, it will be assumed that the loss and smoothness coef- 
ficients are constant over the entire surface; this means especially 
that they are independent of the angle of incidence. Otherwise, numer- 
ical integration is called for. 


W. Heybey, "On Radiation Pressure and Its Effects on Satellite Motion," 
NASA TM X-53462, April 1966. 
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III. THE SPHERE 


If the sphere's center is at the origin, the parameters a and t 
can be chosen as the angles 9 and a in the system of spherical coordin- 
ates (figure 1). The representation (11) assumes the form 


x = R sin a cos cp 
y = R sin -3 sin cp 
z = R cos -a 


where R is the radius of the sphere. All points of the sphere are 
accounted for if the parameters cp and -3 move in the intervals 


0 g 9 g 2it 

0 s ^ it. 

The coordinate line% 9 = const, and -a * const, are the circles of longi-' 
tude and latitude, respectively. 

Application of expressions (12) and (13) yields 


n^ = -sin •a cos 9, n^ = -sin a sin 9, n 3 = -cos a 
N - R 2 sin a. 


Then 


Let the parallel light arrive directly from above: 
s = -k. 

cos a * n • s - cos a. 


d 2 F ■ p Q R 2 sin a da d9 


The differential (10) becomes 

- k (1 - q/c) cos a - 2 k + q cos 2 a 

(± sln ^ cos <p + i sin 4 sin cp + £ cos a)] . (16) 


8 ) 


9 


The upper half of the sphere only is illuminated. While one has 
to integrate over the full natural interval 0 g cp g 2jt, that of -3 must 
be curtailed into 0 g 3 g jt/2, because a shadow boundary exists at 
-9 = it/ 2. 

On integrating over cp, the _i- and j_- constants of the force vanish, 
so that 


dF = -2n p R 2 sin 3 

— o 


(1 - q k) cos 3 + 2 K 



+ q cos 



cos 


3 


k d3. 


Integration with respect to 3 yields 

F = -it Pq R 2 (l + 2 k 1 ~ k = -F k . (17) 

This force has the magnitude F and is in the direction -k of the incoming 
light. If m is the mass of the sphere, the center of gravity moves with 
the speed v = (F/m)t in the direction ^ (v = 0 for t = 0). 

For symmetry reasons the line of attack coincides with the zenith- 
nadir-line. Any sphere whose mass center is on that line experiences 
no torque; a homogeneous sphere therefore cannot be set into rotatory 
motion at all through the force exerted by the light, no matter from 
which direction it arrives. 

Suppose, however, that the sphere is divided by a plane, y = y*, 
parallel to the (z,x)-plane such that left of it the mass density is 
smaller than right of it (figure 3). The mass center, G, will then be 
at a distance y > 0 on the y-axis . At t = 0, the sphere is subject to 
a moment 


M = (-j.y g ) x (-Fk) = iy g F, 

which seeks to rotate the body about an ^-axis parallel through the mass 
center and thus moves the geometric center, C, of the sphere through, say, 
an angle Ai|r from its original seat atj§. The force (17) retains magni- 
tude and direction; however, its arm GC is now 



cos AJr - 



sin A\|f, 
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Figure 3. Inhomogeneous Sphere 
so that the moment becomes 


M = iy F 
- —g 


cos A\|r . 


It is still in the ^-direction, which is also the direction of a prin- 
cipal axis of inertia through the mass center G. The motion will there- 
fore be a mere rotation about this axis (no tumbling occurring). Let T 
denote the moment of inertia relative to the axis of rotation. 

With a finite angle \|r, that is, at a given time t during the motion, 
the absolute value of the moment, M = y. F cos \|r, appears in the equa- 
tion of motion 


dfjfc _ 


A cos \|r 


where A = ygF/T is a positive quantity. Putting d\|r/dt « w, this equa- 
tion may be written as 
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it ^ = 2a ^ cos ♦ 


so that a first integral, with \|r = 0 and go = 0 for t = 0, becomes 


co 


± n/2a" 


sin \[r 


dt ’ 


Further integration shows t to be an elliptic integral similar to the 
integral that occurs with the motion of a pendulum. In writing the 
solution down let it be recalled that the Jacobian elliptic function 
sinus amplitudinis (sn z) has the real period 4K and takes on the 
values 0,1, 0,-1 for z = 0,K,2K,3K, respectively. With its aid the 
dependence of the angle \|r on time may be given as 


sn 2 (n/a t) 

sin \|r = . 

2 - sn 2 (n/a t) 


The angular velocity then becomes 


+ = 


co 



2A 

sn 2 (\]~A t) 


sn 


<S/A t). 


The modulus characteristic for the function here is k = 1 /n/T, so that 


K = 1.8541, 


as can be found from a table of the complete elliptic integrals. 

Now, at t = 0, \Jr = 0 and co = 0. As time proceeds, the angular 
velocity attains positive values up to a maximum \T2A at \[r = 90° 

(4T t = K), when the sphere 1 s original right-hand extremity is now on 
top and the torque M has become zero (the mass center being on the line 
of attack; cos \|r = 0) . Inertia carries the motion on; a torque then 
appears in the negative _i-direction and causes deceleration. The 
angular velocity decreases to the value co = 0 at \|r - 180° (n/X t = 2K) , 


12 


where the "negative" torque has its greatest absolute value and succeeds 
in reverting the sign of the angular velocity. A return half-cycle 
follows which takes the velocity through negative values, until at 
vA t = 4K the initial position and velocity are attained a second time. 

Since the quantity 

A = • Jtp Q R 2 ^1 + 2/c ^ 

can be expected to be very small, especially with appreciable sphere 
mass, this swinging motion will often be slow. If, for example, n/a 
is of the order 10 -5 1/sec, the time needed to accomplish a full cycle 
would be 4K x 10 s sec or somewhat more than 8 days. 

In addition, other small torques may be operative (through micro- 
meteoric impact, the gravity gradient, aerodynamic forces) which may 
lead to a slow gyroscopic (tumbling) motion after all. 

If the force has the general direction s^ (instead of -k) , we may 
introduce the k-direction in such a way that £ is in the (y,z)-plane. 

It is then not hard to show that qualitatively the motion is the same 
as before (assuming again that \|r = w = 0 for t = 0) . In the exceptional 
cases where s^ = ± j, no rotatory motion at all is initiated for lack of 
a moment at t = 0. 

Expression (17) shows that the force takes on a particularly simple 
form for a totally absorbing surface (k = 0). This is to be expected, 
since the sphere then moves under the sole impact of incident light. 

More remarkable is the fact that the same result obtains if the sur- 
face is a perfect mirror (q = 1, K arbitrary). Quite regardless of 
what part of the energy goes into the reflected light, the force is 
always due to incidence alone. One can indeed see from expression (16) 
that the vertical (k-) components of the forces d 2 F r cancel at surface 
points that lie 45° apart on a meridian. The lateral components cancel 
at points 180° apart on a circle of latitude. The net force exerted 
through specular reflection is therefore zero. 
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IV. THE CIRCULAR CONE 


In figure 4 it is assumed that the light, arriving from the left 
and back, determines the (x,y)-plane with the cone (y-) axis. Its 

2, JL 



direction is then given by the unit vector 


s = i sin p + j_ cos p 


where p is counted a positive angle (negative values would assign a 
backward component to the light direction with no essential consequences 
for the description). Let us consider the case of full illumination: 


P ^ to. 


The parametric representation of the curved surface may be written as 


x = Ra sin t 
y = La 
z = Ra cos t 
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where 


0 g a S 1, 0 £ x £ 2it 


are the natural intervals needed for full coverage of the surface. The 
coordinate lines a = const., t = const, are circles normal to the cone 
axis and the generating lines, respectively. 

Since 


|y z I 


Z X 


x y 

^ Q 

N 

H Q 

= -LRa sin x. 

a a 

Z X 
T T 

= R 2 a, 

a J a 
x y 

T J T 


= -LRa cos 


T, 


it follows from the general expressions (13), (12) and (15) that 


N = Rcr n/l 2 + R 2 = RLg 

cos to 

nj. = -sin t cos to, n 2 = sin to, n 3 = -cos t cos to 

cos a = -sin t cos to sin p + sin to cos p. 

With these results and the area differential (14) the force element 
(10) becomes 

p 0 R LcJ f 

d 2 F = ' CQ ~ co ~' da dr -j (1 - qk) (_i sin P + j. cos p)(-sin x cos to sin p 

+ sin u cos p) + 2/c(-_i sin x cos to + j_sinto-k cos x cos to) 

• - - + q (-sin x cos u sin p + sin to cos p) 2 
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In integrating over the natural a- interval the factor 1/2 appears. 
Since the parameter t also moves in its full interval 0 S t S 2jt (no 
shadow boundary being admitted), all r-quadratures leading to trigono- 
metric functions alone result in zero. The final outcome is less 
involved than might have been expected; 


F = xp R 2 -j_i sin p cos (3 (1 + qk cos 2w) 


+ i 


cos 2 p + — ^ ^ K - Kq (sin 2 f3 sin 2 u + 2 cos 2 |3 cos 2 co)J 


(18) 


The cone length L appearing in the force differential has been replaced 
here by R cotg oj. 

It is seen that with q = 1, expression (18) still contains terms 
involving the loss coefficient K and that therefore, other than with 
the sphere, the force continues to depend on the quality of the reflec- 
tion. It is not equal to the force through incidence alone; 


F i = 3tP Q R 2 cos (3 (jL sin |3 + cos p) (19) 

except in the trivial case where there is no reflection (K = 0). Its 
absolute value is often smaller than that of f^, especially with slender 
cones, with which the force through specular reflection has a strong 
component opposite to the axial component of the incident force. Con- 
sider zero incidence where 



R 2 


while, by expression (18) 


F| = Jtp Q R 2 


1 + K 


2 + q _ 


2q cos^co 


■)]• 


When 


COS 2 OJ > 


2 + 
6q 


a 


F | < F-l . In the ideal case q = 1 the condition becomes w < 45°. 
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Clearly, reflection not always strengthens the force as one might 
be led to believe by incorrectly generalizing flat plate results. 


Fairly extended calculation is necessary to locate the center of 
pressure. It is found at 


1 + 




x* = z* = 0, y* = L 


1 + q/c cos 2oo 


Its position does not depend on the angle p, but it is variable with 
the cone angle. When the cone length is fixed, it moves away from 
the tip with increasing base radius; if go reaches a value « 35°15' 
(cos 2go = 1/3), it resides at base center. 

For a homogeneous cone the mass center is at 



The stability condition, y* > y g , tells that most cones are stable 
except the slender ones for which 


1 + 


a* 

3 


1 + 2/C cos 2io 



( 20 ) 


As long as q/c < 0.6, the cone motion is stable, in particular with a 
non-reflecting surface (/c = 0) and with completely diffuse reflection 
(q = 0). Ideal specular reflection (/c = q = 1) renders a cone unstable 
provided its tip angle is sufficiently small (sin go < 1/3, go < 19.5°). 

The torque about the center of pressure is found from 


M = - j (y - y*) x F 

o 


which expression gives 


M = itp Q R 2 L sin p cos 



(1 + q/c cos 


2 go) 



( 21 ) 
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This formula, as the force expression (18), is valid only with a fully 
illuminated cone, i.e., for p I u. Once the angle p is inside this 
region, it will remain there provided the motion is stable. Otherwise, 

P will eventually exceed co. The much more complicated expressions which 
arise with the presence of a shadow boundary on the cone simplify for 
a non-reflecting surface; one finds that, with p & oo, the motion is 
stable here as long as tan w cotg p > 0.17. 

In the case of the Apollo capsule where co « 33°, this complication 
is not likely to occur*, although it is not a homogeneous body. Once 
the tip points in the general direction of the sun so that the conical 
part of the surface is entirely exposed to sunlight, the torque pro- 
vided by the latter should tend to align the cone axis with light direc- 
tion. It would be expected that the same tendency prevails if a narrow 
triangular shadow zone exists on the surface. However, with a wider 
such zone, a detailed mathematical investigation is called for, which 
also would have to allow for the forces that attack the illuminated 
parts of the base cap. 

In order to study the behavior of the light force acting on a fully 
illuminated cone„a computer program** was run which assumed lossless 
reflection ( K =1). The angle co was taken as 33° (Apollo capsule). 

The remaining parameters then are the smoothness coefficient q and the 
angle p (^ to) of overall incidence. The results can be summarized as 
follows : 

(a) p = const. The force increases with the smoothness 
coefficient decreasing from 1 to zero (in the approximate 
ratios 1:3, 3:5 at p = 0°, p= 33°). Its inclination 
towards the cone axis decreases (except at p = 0 where it 
is always zero). The maximum decrease occurs at p = 33° 
(from 45.75° to 18.44°). In short, increasing diffusivity 
renders the force stronger and less inclined towards the 
cone axis . 

(b) q = const. The inclination of the force increases with 

P (from 0 to 18.44° at q = 0, from 0 to 45.75° at q = 1). 
Its magnitude decreases if q <0.5 (in the ratio 7:6 for 
q = 0); it increases if q >0.5 (in the ratio 2:3 for 
q = 1). With q = 0.5 there is an insignificant force 
maximum at p « 25°, the whole variation amounting to less 
than 1 percent in the full range 0 ^ p ^ to. 

“ft 

There is no value of q K 1) that would satisfy condition (20), as 
cos 2 (a) » 0.4. 




Set up by James E. Mabry. 
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With the opposite extreme of K (k = 0) , the behavior of the force 
does not depend on co and is directly seen from expression (19). Its 
angle with the cone axis is always 3 , and its absolute value decreases 
as 3 moves from 0° to co. 

If 3 exceeds to, expression (19) for to = 0 must be replaced by 


F = J\ = P q R 2 (i sin 3 + j. cos 3) 


r it + 2 t 


— cos 3 + sin 3 cotg to cos T g 


where the relation 


sin t = tg to cotg 3 
s 

defines the two straight shadow boundaries t = t s and t = it - 2x g . 

In the limiting cases of the half and fully illuminated cone 
(3 = it/ 2 , 3 = w ), one has 


ZlfS-K /2 = p o R2 cotg “ 

I F I = itp R 2 cos to. 

It is seen that if sin to < 1/ it the force at 3 = it/ 2 is stronger than 
that at 3 = u; that is, with sufficiently slender cones the decrease 
of the force with increasing 3 will not continue indefinitely up to 
3 = it/2. It can be shown that in such circumstances, the force goes 
through a minimum followed by a maximum and ends up decreasing at 
3 = it/2. In the special case where the force magnitude here is iden- 
tical with that at p = a) (sin co = 1/it, co 17.5°), the minimum was found 
at 3 ** 19.5°, the maximum at 3 « 58°. 

With co = 33°, the magnitude of the force on a non-reflecting cone 
decreases through the entire 3 - interval from 3 = 0 to 3 = it/2. From 
unity at 3 “ 0 it: is down to 0.84 at 3 = co, and to 0.49 at 3 ■ rt/2. 
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V. CONCLUDING REMARKS 


Diffuse re-emission was thought in the foregoing to be governed by 
Lambert* s cosine law. This law furnished the model from which to cal- 
culate the overall diffusion force acting on an elemental surface. From 
the intensities re-emitted in the several directions, corresponding light 
forces were determined and added up vectorially. 

Other approaches to the problem of diffuse reflectivity can be 
imagined . 

As an example, one may assume that an infinitesimal half-sphere 
erected over the elemental surface is uniformly illuminated by the 
diffusely reflected light. Contrary to Lambert* s rule, its intens- 
ity then would be the same in all directions. This entails that 
a greater part than before of the re-emitted energy is caught up in 
cancelling **s tresses** parallel to the elemental surface. Computation 
shows that the factor 2/3 in expressions (4) and (9) is to be replaced 
by 1/2, so that instead of (1 - q)/3 in the fundamental expression (10) 
one has to write (1 - q)/4. 

In both these approaches fictitious model forces are treated as if 
they were actual forces; the factors 2/3 and 1/2 are uncertain to this 
extent. From the random distribution of the surface irregularities one 
would conclude that the integrated force actually exerted on the sur- 
face element in diffuse reflection should be essentially in the direc- 
tion of the interior surface normal (as in fact.it turned out to be 
with the two models). Without going into any more detail, one might 
argue that the amount of energy not reflected specularly is reflected 
diffusely, thus converting the factor 2/3 in expression (4) into unity 
and the term (1 - q)/3 into (1 - q)/2. 

From all this it would seem that the **dif fus ion** portion of the 
force differential (10) is likely to be proportional to a factor between 
(1 - q)/2 and (1 - q)/4. 

To illustrate the differences involved the forces arising with 
K - 1 on the fully illuminated cone with half-angle go = 33° were com- 
puted using the largest factor (1 - q)/2 in order to compare them with 
those previously found for (1 - q)/3. In the worst case (q = 0; exclus- 
ively diffuse reflection) the force magnitudes are multiplied by nearly 
the same factor 1.2 in the entire range 0 ^ p ^ oo, whereas the force 
inclinations toward the cone axis decrease by a rather constant factor 
of again ^1.2. A tendency of this kind is in keeping with what an 
earlier result should suggest to expect from increasing diffusivity; 
the force on the cone becomes stronger and less inclined. Its moment 
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about the origin, however, and thus the center of pressure, are inde- 
pendent of speculations on the magnitude of the diffusion part, because 
with full illumination the latter's line of attack is the cone axis. 

The diffusion torque therefore is zero irrespective of the choice of 
n in (1 - q)/n and, by the way, of all other characteristic quantities 
such R, L, q, K. 
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